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The  AFOSR  program  at  UCSD  contains  mainly  two  areas  of  optical  information  processing 


research: 


(i)  Optical-analog/electronic  hybrid  computing 


On  hybrid  optical-analog/electronic  computing,  we  have  performed  several  new  studies  for 
optical  image  processing,  optical  pattern  recognition  and  optical  solution  of  partial  differential 
equations.  For  example,  for  image  processing  we  designed  and  fabricated  space  variant  filters  to 
optically  implement  numerous  space-variant  transformations  (e.g..  Hough  transform  for  detection 
of  high-order  parametric  curves,  coordinate  transforms  to  detect  rotation  and  scale  invariant 
features  of  images,  etc.).  For  spike  removal  from  noisy  images,  we  have  developed  a  new  parallel 
algorithm  suitable  for  optical-analog/electronic  hybrid  implementation.  This  parallel  algorithm  is 
space-variant  and  performs  better  than  the  space-invariant  low-pass  filter  and  the  time  consuming 
median  filter.  We  also  studied  a  single  beam  interferometer,  which  effectively  subtracts  an 
exponentially  weighted  history  of  the  input  from  the  current  value,  thus  functioning  as  a  time- 
differentiating  (or  novelty)  filter.  The  single  beam  interferometer  uses  signal  beam  depletion  from 
photorefractive  amplification  of  noise  (fanout)  in  a  specially  designed  and  cut  BaTiO  3 
photorefractive  crystal  (with  input  faces  cut  parallel  to  (100),  (Oil)  and  (Oil)  crystallographic 
planes). 

For  optical  pattern  recognition  we  have  completed  the  comparison  study  of  statistical  pattern 
recognition  algorithms  by  investigating  the  performances  of  their  discriminant  filters  when 
implemented  by  phase  only  (PO)  filters.  We  found  that  for  each  of  the  pattern  recognition 
algorithms  studied  (e.g.,  least-square  linear  mapping  technique  (LSLMT)  simplified  LSLMT 
(SLSLMT),  Hotelling  Trace  Criterion  (HTC),  Foley-Sammon  (FS)  transform,  and  Fukunaga- 
Koontz  (FK)  transform)  the  peak  of  the  correlation  intensity  increases  and  the  size  of  the  correlation 
spot  decreases  when  PO  filter  is  used.  The  signal-to-noise  ratio  wa«  improved  for  the  LSLMT, 
SLSLMT  and  the  FK  filters,  but  deteriorated  for  the  HTC  and  the  FS  filters.  Continuing  our  search 
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for  new,  more  robust  optical  pattern  recogni'ion  algorithms,  we  have  formulated  a  nonlinear 
transformation  based  on  minimization  of  a  special  energy  function,  (i.e.,  self-organization  neural 
network).  In  contrast  to  those  supervised  learning  algorithms  based  on  data  of  known  statistics,  this 


unsupervised  learning  algorithm  does  not  require  a  priori  knowledge  of  the  statistics  of  the  input 
data.  To  implement  the  nonlinear  transformation  for  optical  pattern  recognition  we  are  constructing 
a  hybrid  optical/electronic  system  based  on  random-phase  coded  matrix-tensor  multiplier,  which 
can  provide  in  parallel  and  in  real  time  a  large  number  (e.g.  —  1 06  —  1 09 )  of  adaptive  optical 
connections  and  therefore  allow  operation  on  large  size  images. 

For  the  optical  solution  of  partial  differential  equations  by  a  confocal  Fabry-Perot  system  with 
optical  feedback  and  photorefractive  gain,  we  have  shown  that  the  analog  accuracy  of  the  optical 
solutions  to  partial  differential  equations  can  be  improved  by  implementing  an  iterative  error 
correcting  algorithm. 

(ii)  Digital  optical  computing  and  nonlinear  optical  devices: 

On  digital  optical  computing  research,  we  have  been  developing  nonlinear  optical  devices  and 
studying  their  uses  in  parallel  architectures  for  implementation  of  parallel  algorithms.  The 
nonlinear  optical  devices  are  based  on  integrating  silicon  with  PLZT.  Using  LPCVD  techniques  to 
deposit  polycrystalline  Si  onto  the  electro-optic  PLZT  substrate  and  then  recrystalizing  the 
polysilicon  using  Ar+  laser  to  produce  single  crystal  silicon  grains,  we  have  fabricated  and  tested  a 
12  x  12  electrically  addressed  spatial  light  modulator  (SLM)  array  and  a  16  x  16  optically 
addressed  SLM  array. 

To  study  optical  parallel  computing  architectures,  we  investigated  different  interconnection 
topologies  and  have  generated  holographic  optical  elements  to  implement  fully  interconnected  and 
hypercube  interconnected  topologies. 


1.  Introduction 
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The  AFOSR  program  at  UCSD  consists  of  two  areas  of  optical  information  processing:  (i) 
optical-analog/electronic  hybrid  computing  for  optical  image  processing,  optical  pattern  recognition 
and  solution  of  partial  differential  equations,  and  (ii)  nonlinear  optical  devices  and  digital  optical 
computing.  The  optical-analog/electronic  digital  hybrid  system  consists  of  an  analog  optical 
processor  with  input/output  interfaces  to  a  microcomputer.  The  analog  optical  processor  is 
employed  to  perform  time  consuming  computations,  while  the  logical  decisions  and  controls  are 
provided  by  an  electronic  microcomputer.  For  image  processing  we  have  implemented,  with  hybrid 
systems,  numerous  space-variant  transformations  (e.g.,  Hough  transform  for  detection  of  high-order 
parametric  curves,  coordinate  transform  for  rotation  and  scale  invariant  feature  extraction,  etc.).  For 
spike  removal  from  noisy  images  we  have  developed  a  new  parallel  algorithm  for  optical 
implementation,  whose  performance  is  superior  to  that  of  low-pass  or  median  filters.  For  motion 
detection,  image  clutter  removal  and  novelty  filtering  we  have  constructed  and  experimentally 
evaluated  a  time-differentiating  (or  novelty)  filter,  which  consists  of  a  single  beam  interferometer 
based  on  a  specially  designed  BaTiO  3  photorefractive  crystal  (with  input  faces  cut  parallel  to  (100), 
(Oil)  and  (Oil)  crystallographic  planes).  For  optical  statistical  pattern  recognition  we  studied  the 
implementation  of  discriminant  functions  of  our  phase-coded  processor  using  phase  only  filters. 
We  also  formulated  a  new  nonlinear  transform  for  unsupervised  optical  partem  recognition,  which 
is  presently  being  implemented  by  a  hybrid  optical/electronic  system  based  on  random-phase  coded 
matrix -tensor  multiplier.  Finally,  for  the  optical  solution  of  partial  differential  equations  with  the 
hybrid  system  based  on  confocal  Fabry-Perot  processor,  we  have  shown  that  the  analog  accuracy  of 
the  optical  solutions  can  be  improved  by  implementing  an  iterative  error  correcting  algorithm. 

On  the  digital  optical  computing  research,  we  have  been  developing  nonlinear  optical  devices 
and  studying  their  uses  in  parallel  architectures  for  implementations  parallel  algorithms.  Nonlinear 
optical  devices  based  on  integrating  silicon  with  PLZT  were  fabricated  and  tested  (e.g.,  a  12  x  12 
eiectricaiiy  addressed  spatial  light  modulator  and  a  16  x  16  optically  addressed  spatial  light 
modulator).  In  our  attempt  to  construct  digital  optical  processing  systems  we  have  fabricated 


optical  holograms  for  implementation  of  architectures  of  fully  and  hypercube  interconnection 
topology. 
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2.  Hybrid  (Optical-Analog/Eledronic  Micro-)  Computer 


2.1.  Optical  Image  Processing 


2.1.1.  Optical  space-variant  filters  for  image  transformations 

In  the  study  of  digital  image  processing  algorithms,  which  would  benefit  from  optical 
implementation,  we  have  constructed  a  computerized  optical  system  for  design  and  generation  of 
space-variant  holographic  filters  (see  Fig.  1)  [1].  Employing  this  system,  space-variant  filters 
containing  256  x  256  sub-holograms  have  been  generated  in  649F  and  dichromatic  gelatin  plates 
for  optical  Hough  transforms  of  high-order  parametric  curves  and  for  coordinate  transforms.  These 
filters  have  been  used  in  the  real  time  space-variant  optical/electronic  processor  of  Fig.  2  to  detect 
the  three  parameters  of  a  circle  and  the  four  parameters  of  an  ellipse;  the  experimental  performance 
is  summarized  in  Fig.  3  and  Fig.  4  respectively.  Furthermore,  to  perform  rotation  and  scale 
invariant  image  transformation,  we  have  generated  a  in  r-G  coordinate  transform  filter  for  images 
of  size  256  x  256  pixels.  The  experimental  results  on  performance  of  this  filter  are  shown  in  Fig.  5 


and  Fig.  6.  Details  can  be  found  in  Ref.  1. 


2.1.2.  Parallel  optical  algorithm  for  image  spike  noise  removal 

We  have  developed  a  new  parallel  algorithm  for  image  spike  noise  removal  suitable  for 
optical  processing  implementation.  The  new  algorithm  is  space-variant  and  performs  better  than 
the  simple  space-invariant  low-pass  filter  used  to  reduce  noise  in  low-level  image  processing,  (see 
Figure  7(b)),  and  the  time-consuming  median  filter  (see  Figure  7(c)).  The  space-variant  algorithm 
for  spike  noise  removal  can  be  applied  to  a  large  size  image  and  is  based  on  comparing  the  intensity 
of  a  pixel  /(/,_/)  with  the  average  intensity  of  its  nearest  neighbors  (/(/,;)).  If  \I(iyj)-I(i,j)\  is 
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greater  than  the  local  variance  of  the  nearest  neighbors,  then  we  replace  I(i.j)  by  I(i,j),  otherwise 
the  original  value  of  / (i,j)  is  preserved.  Computer  simulation  results  (see  Figure  7(d))  show  that 
the  proposed  algorithm  removes  spike  noise  without  blurring  the  image.  Table  I  summarizes 
quantitatively  the  perfonnance  of  the  proposed  algorithm  and  compares  it  with  the  conventionally 
used  spike  noise  removal  algorithms  (i.e..  median  and  low  pass  filter). 


2.1.3.  Optical  time-differentiating  filter  for  image  motion  analvsis 


We  have  demonstrated  a  single  beam  interferometer  which  effectively  subtracts  an 
exponentially  weighted  history  of  the  input  from  the  current  value,  thus  functioning  as  a  novelty 
filter  [2],  The  filter  uses  signal  depletion  from  photorefractive  amplification  of  noise  (fanout). 
Input  which  is  stationary  for  longer  than  the  photorefractive  response  time  is  strongly  depleted,  but 
any  portions  of  the  input  which  change  their  phase,  amplitude,  wavelength  or  polarization  are 
immediately  transmitted.  The  gain  material  used  was  a  crystal  of  photorefractive  barium  titanate 
with  its  surfaces  on  the  (1,0,0),  (0,1,1)  and  (0,1,1)  crystallographic  planes,  so  that  the  optic  axis  lies 
at  45  degrees  to  the  input  face,  allowing  access  to  the  peak  gain  with  normal  signal  incidence  [3]. 
This  crystal  provided  higher  gain  and  faster  response  than  the  normally  cut  material.  To  examine 
the  time-differentiating  filter’s  real  time  operation,  we  used  a  hybrid  optical/electronic  system  with 
input/output  interfaces:  a  liquid  crystal  light  valve  is  used  to  convert  a  video  image  into  a  phase 
and/or  amplitude  modulated  input  signal  and  a  CCD  camera  is  used  as  an  optical  to  electronic 
convertor.  The  experimental  results  are  shown  in  Fig.  8.  A  contrast  ratio  of  100  was  observed 
between  moving  and  stationary  images,  with  a  resolution  limited  only  by  that  of  the  light  valve. 
Potential  applications  of  this  filter  include  image  clutter  removal,  motion  detection  and  analysis, 
tracking,  edge  enhancement,  and  image  time  differentiation. 
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2.2.  Optical  Pattern  Recognition 

2.2.1.  Comparative  Study  of  Phase  Only  Filter  for  Statistical  Pattern  Recognition 

We  have  recently  completed  a  study  of  phase-only-filtering  applied  to  discriminant  functions 
provided  by  various  statistical  pattern  recognition  algorithms.  Phase-only  (PO)  filtering  employs 
only  the  phase  of  the  filtering  function  while  the  amplitude  is  set  equal  to  1  at  all  points  on  the 
filter.  We  performed  computer  simulations  in  order  to  evaluate  the  effect  of  PO  filtering  as 
compared  to  amplitude  and  phase  (AP)  filtering.  The  performance  of  PO  and  AP  filtering  were 
evaluated  using  the  following  characteristics:  peak  correlation  intensity,  correlation  spot  size, 
signal-to-noise  ratio  and  optical  efficiency.  We  were  interested  in  measuring  these  characteristics 
for  the  discriminant  functions  provided  by  various  statistical  pattern  recognition  algorithms  that 
have  been  under  investigation  at  UCSD.  The  five  different  algorithms  we  evaluated  in  conjunction 
with  the  PO  filtering  study  were:  the  Least  Squares  Linear  Mapping  Technique  (LSLMT)  [4],  the 
Simplified  Least  Squares  Linear  Mapping  Technique  (SLSLMT)  [5],  the  Hotelling  Trace  Criterion 
(HTC)  [b],  the  Foley-Sammon  Transform  (FS)  [7]  and  the  Fukunaga  Koontz  transform  (FK)  [8]. 
The  training  and  test  data  we  used  consisted  of  grey  tone  images  of  64  x  64  pixels  with  one  byte  of 
grey  tone  data  per  pixel.  Our  training  set  was  composed  of  5  such  images  of  tanks  and  5  of  trucks. 
Our  test  data  set  was  composed  of  4  tank  and  4  truck  images. 

For  each  algorithm  we  created  two  filters  (i.e.,  PO  and  AP)  using  the  two  classes  of  training 
images.  The  test  images  were  then  used  to  compute  the  various  evaluation  characteristics.  The 
simulation  results  are  summmarized  in  Tables  2-6  for  evaluation  of  LSLMT,  SLSLMT,  HTC,  FS 
and  FK  transforms.  The  results  obtained  indicate  that  for  all  algorithms  the  peak  correlation 
intensity  increases  and  the  correlation  spot  size  decreases  when  a  PO  filter  is  used.  The  signal-to- 
noise  ratio  improved  for  the  LSLMT,  SLSLMT  and  the  FK  filters  but  deteriorated  for  the  HTC  ard 
the  FS  filters.  The  optical  efficiency  did  not  show  a  significant  change  in  either  direction. 
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2.2.2.  Nonlinear  Transformations  for  Pattern  Recognition 

The  feature  extraction  process  of  a  partem  recognition  system  can  be  considered  as  a  function 
/  that  maps  a  set  of  input  patterns  x  onto  a  set  of  features  fix)-  The  criteria  usually  used  for 
designing  the  feature  extraction  function  is  to  minimize  an  error  energy  function  £,  which  is 
defined  as  ||/  (.v)-  T(.v||2 ,  where  T(.v)  is  the  desired  output  feature. 

Conventional  pattern  recognition  systems  [4-9]  use  pre-designed  and  fixed  linear  filters  /, 
which  are  implemented  based  on  a  priori  knowledge  of  the  input  patterns.  If  an  adaptive  pattern 
recognition  system  is  desired,  the  function  /  is  implemented  with  a  recursive  algorithm  [10-12], 
where  an  on-line  learning  rule  is  enforced  to  adjust  /  successively.  Usually  the  function  T(.v)  is 
determined  by  the  statistics  of  the  input  pattern.  Therefore,  for  data  of  known  probability  density 
function,  TU)  can  be  considered  given  and  the  adaptive  process  can  be  seen  as  supervised  learning. 
However,  those  systems  will  not  work  well  in  a  changing  environment,  where  the  statistics  of  the 
input  patterns  are  usually  not  predictable,  or  do  not  obey  the  simplified  gaussian  assumption. 
Furthermore,  for  some  complicated  patterns,  due  to  nonlinear  characteristics,  patterns  may  not  be 
well  classified  by  the  most  frequently  used  linear  discriminant  function  in  the  feature  domain. 
Therefore,  nonlinear  transformations  or  nonlinear  feature  extraction  processes  should  be  employed 
to  extract  more  effective  features  in  order  to  gain  higher  accuracy  and  reliability  in  pattern 
classification. 

We  have  already  shown  that  an  adaptive  pattern  recognition  system  may  be  put  in  terms  of 
neural  network  systems  as  an  adaptive  hetero-associative  memory  [12].  The  learning  process  of  a 
hetero-associative  memory  is  equivalent  to  the  supervised  on-line  training  of  an  adaptive  pattern 
recognition  system,  while  recalling  information  from  a  hetero-associative  memory  is  equivalent  to 
the  feature  extraction  process  of  an  adaptive  pattern  recognition  system.  Both  the  adaption  rule  of  a 
pattern  recognition  system  and  the  learning  rule  of  a  hetero-associative  memory  are  methods  that 
minimize  the  errors  between  the  actual  outputs  and  the  desired  outputs. 

To  overcome  the  problems  of  conventional  pattern  recognition  systems  as  stated  above,  we 
have  been  developing  a  neural  network  which  utilizes  the  concept  of  self-organization  [11,13,14]  to 


optimize  the  solution  without  having  input  statistics.  The  proposed  neural  network  is  based  on  the 
minimization  of  a  special  energy  function,  £'  =  ||/2(/i(-0)- v||2,  where  both  fy  and  fi  are,  in 
general,  nonlinear.  In  fact,  the  energy  function  is  a  mean-square  difference  (error)  between  the  true 
pattern  ,v  and  the  estimate  pattern  v  =  fi(f  i(.v)).  Without  knowing  any  input  statistics,  the  network 
will  self-organize  its  interconnections  by  adaptively  modifying  their  strength  in  a  direction  that 
minimizes  the  energy  function  £' .  This  energy  minimization  process  will  result  in  optimum 
transformations  fy  and  /;,  so  that  the  resultant  self-organizing  feature  vector  fy  (a  )  will  optimally 
describe  the  input  data  statistics.  This  is  different  from  the  self-organization  algorithm  developed 
from  the  gradient  method  [11,13]  and  the  algorithm  based  on  cooperative-competitive  learning 
[14]. 

The  energy  minimization  pracess  for  self-organization  will  require  the  development  of 
techniques  capable  of  finding  the  desired  optimum  nonlinear  transformations  fy  and  /? ■  In 
supervised  learning  one  of  the  widely  used  techniques  is  the  delta  rule  [9].  We  are  studying  the 
extension  of  the  delta  rule  as  applied  to  the  minimization  of  the  energy  function  E'  for  unsupervised 
learning.  Another  technique  used  in  unsupervised  learning  is  the  simulated  annealing  algorithm 
which  is,  in  fact,  a  stochastic  approach.  The  simulated  annealing  algorithm  is  well  suited  for  the 
minimization  of  the  proposed  energy  function,  and  therefore  is  also  under  investigation. 

To  implement  the  discussed  algorithm  optically  we  are  presently  constructing  an  optical 
neural  network  architecture  based  on  a  real  time  random  phase  coded  matrix-tensor  multiplier 
interfaced  to  a  microcomputer.  This  architecture  (see  Fig.  9)  ■will  be  capable  of  storing  a  very  large 
number  of  optical  connections  (~106  -  109)  using  a  photore tractive  crystal.  All  these  connections 
can  be  accessed  in  parallel  and  in  real  time  to  provide  the  output  matrix  for  electronic  thresholding 
and  calculation  of  a  correction  term  for  the  tensor  according  to  the  learning  rule.  The  correction 
term  will  be  used  to  update  the  optical  interconnection  stored  in  the  photorefractive  memory 


material. 
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2.3.  Solution  of  Partial  Differential  Equations 

The  Bimodai  Error  Reduction  (BER)  algorithm  is  used  for  improving  the  accuracy  of  the 
solutions  obtained  from  a  moderate-accuracy  processor  [15.16].  Our  investigation  involves  the 
experimental  application  of  this  algorithm  to  satisfy  the  desired  accuracy  requirements  on  an 
electronic/optical  hybrid  processor  being  constructed.  The  optical  portion  of  this  hybrid  processor 
performs  the  analog  solution  of  Poisson's  equation  for  a  two-dimensional  input  image.  The 
electronic  portion  updates  the  input  image  and  performs  support  functions.  We  have  theoretically 
analyzed  the  error  correction  algorithm  applied  to  Poisson's  operator,  which  led  to  an  expression 
for  the  condition  of  convergence.  Also,  a  computer  simulation  of  the  hybrid  (electronic/optical) 
implementation  of  the  algorithm  was  performed  to  confirm  the  theoretical  analysis.  The  theoretical 
development  considered  additive  error  on  the  operator  in  the  analog  processor  and  on  the  input 
image  presented  to  the  analog  processor.  Both  error  functions  were  assumed  to  be  spatially  and 
temporally  constant  in  this  development.  Conditions  for  convergence  of  the  algorithm  under  these 
assumptions  were  determined  and  found  to  include  the  expected  operating  conditions  of  the  hybrid 
(analog/digital)  processor  being  constructed. 

A  simulation  of  the  BER  algorithm  as  implemented  on  the  hybrid  processor  confirmed  the 
theoretical  results  (see  Table  7).  Refinement  of  the  solution  to  8-bit  accuracy  is  obtained  in  four 
iterations.  Next  we  enhanced  the  computer  simulation  to  more  realistically  represent  the 
cha  uzteristics  of  the  hybrid  processor.  For  example.  Table  8  shows  the  results  obtained  when 
using  a  sinusoidal  (i.e.,  multi-grey  level)  input  image  and  a  sinusoidal  error  function  for  the 
operator.  Again  convergence  of  the  algorithm  is  obtained.  In  particular,  we  see  achievement  of 
better  than  8-bit  accuracy  in  the  final  solution  after  only  four  error  correcting  iterations. 


3.  Nonlinear  Optical  Devices  and  Digital  Optical  Cnmputin 


3.1.  SiTLZT  Opto-Electronic  Devices  for  POP 

Using  well-established  thin  film  deposition  techniques  to  deposit  polycrystalline  Si  onto  the 
electro-optic  PLZT  substrate  and  then  recrystallize  the  polysilicon  using  an  Ar+  laser  to  produce 
single  crystal  silicon  grains,  vve  have  fabricated  a  12  x  12  electrically-addressed  spatial  light 
modulator  (SLM)  array  and  a  16  x  16  optically-addressed  SLM  array  [17J.  The  schematic  diagrams 
of  the  electronic  circuits  for  the  two  devices  are  shown  in  Fig.  10. 

Each  pixel  of  the  ESLM  array  consists  of  a  driving  transistor,  a  reset  line  and  the  E-0 
modulator.  .An  electrical  input  drives  the  NMOS  transistor,  the  output  of  which  in  turn  controls  the 
modulator.  In  experiments  performed  in  the  last  year,  we  have  achieved  excellent  transistor 
performance  characteristics  which  are  shown  in  Fig.  11  [18].  Figure  11a  shows  the  I-V 
characteristics  of  the  NMOS  transistor  and  Fig.  11b  shows  the  time  response  characteristics  of  the 
ESLM.  We  have  achieved  breakdown  voltages  of  48V  which  is  at  least  a  factor  of  2  greater  than 
that  of  previously  reported  thin  film  transistors  [19].  The  transconductance  was  200 fis  which  then 
yielded  a  calculated  value  of  550  cm2/Vs  for  the  mobility.  We  are  at  present  in  the  process  of 
measuring  the  storage  time,  the  leakage  current  and  the  optical  modulation  depth  of  these 
modulator  arrays. 

The  performance  characteristics  of  the  OSLM  array  are  shown  in  Fig.  12.  In  this  device  we 
have  demonstrated  the  fabrication  of  a  CMOS  device  on  an  electro-optic  substrate.  Breakdown 
voltages  of  28v  were  achieved  for  the  NMOSFET.  The  transconductance  was  \00ji s/V  and  the 
mobility  w-as  442  cm“/Vs  [19].  The  major  improvement  of  the  performance  characteristics  resulted 
from  the  optimization  of  the  process  steps  and  the  improved  crystal  grain  growth  obtained  using 
laser  recrystallization. 
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3.2.  Digital  Optical  Processing 

To  study  optical  parallel  computing  architectures,  we  investigated  different  interconnection 
topologies  and  have  generated  holographic  optical  elements  (HOEs)  that  implement  ta)  fully 
interconnected  and  (b)  hypercube  mterconnected  topologies  [20],  These  HOEs  were  space-variant 
holograms  designed  and  generated  by  the  computerized  system  of  Fig.  1,  and  were  recorded  in 
dichromatic  gelatin  to  provide  diffraction  efficiency  of  better  than  80%.  Reconstructions  from  one 
sub-hologram  (PE)  for  fully  connected  and  one  for  a  hypercube  interconnected  architecture  are 
shown  in  Figs.  13  and  14,  respectively. 
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Table  4.  Comparative  Correlation  Results 


Comparative  Correlation  Results 


Table  7. 


Comparison  of  results  for  the  BER  algorithm  obtained  from  a  theoretical  description  of  the 
algorithm  and  from  a  computer  simulation  of  a  hybrid  (electronic/optical)  implementation  of  the 
algorithm.  The  initial  input  was  a  constant  intensity  two-dimensional  image;  relative  operator  error 
(25%)  was  constant  over  the  input  image  domain.  Experimentally  realizable  lesults  will  be  limited 
by  the  accuracy  of  the  electronic  processor. 
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Results  from  a  computer  simulation  of  a  hybrid  (electronic/optical)  implementation  of  the  BER 
algorithm  using  a  sinusoidal  input  image  function  and  a  sinusoidal  relative  operator  error  function 
(max  25%). 
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Schematic  of  the  hybrid  system  for  producing  computer  generated  space-variant 
holographic  filters.  A  matrix  of  256  by  256  holograms  of  different  PSF  is  recorded 
with  this  system.  The  following  focal  length  (f)  in  mm  and  F-number  (F#)  lenses  are 
employed  f  \  =  100,  fi  —  20,  fi  =  360,  ft  =  80  of  F#  =  1 .6,  =  80  of  F#  =  1 .6,  fa  =  80 

of  F#  =  1.6,  /7  =  20,  /g  =  120,  /9  =  120. 
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Real  time  space-variant  optical  processor  (Top  view).  The  input  pattern  is  displayed 
on  the  CRT/LCLV  and  imaged  on  the  space-variant  holographic  filter.  A  camera  is 
placed  in  the  reconstruction  plane  of  the  holograms. 
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Fig.  5.  In  r-0  coordinate  transform  of  a  rotation  variant  pattern,  a)  Original  input  image  of 
an  off-center  circle,  b)  Optical  coordinate  transform  fnr  -6  of  the  input  pattern,  c) 
Optica!  coordinate  transform  of  the  input  pattern  rotated  by  an  angle  of  45°  around  the 
origin  of  the  coordinate  system,  d)  Optical  coordinate  transform  of  the  input  pattern 
rotated  by  an  angle  of  90°  around  the  origin  of  the  coordinate  system. 


fn  r-0  coordinate  transform  of  another  scale  variant  pattern,  a)  Input  image  of  an 
airplane  with  original  scale  1.  b)  Optical  coordinate  transform  of  the  image  in  (a),  c) 
Input  image  of  (a)  scaled  by  a  factor  of  a  =  0.7.  d)  Optical  coordinate  transform  of  die 
image  in  (c).  The  transformed  pattern  of  (b)  is  shifted  upward  by  an  amount  of  in  a. 


The  transmitted  output  signal  just  after  initial  illumination  of  the  stationary  input  scene. 
After  a  few  seconds,  the  output  vanishes  completely.  Notice  the  model  car  in  the  lower 


The  transmitted  output  just  after  the  model  car  was  moved  to  the  left.  The  car’s  initial 
and  current  positions  are  revealed,  along  with  a  fading  trail  of  the  most  recent  posi¬ 
tions.  The  stationary  components  of  the  scene  remain  dark. 
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Fig.  9.  (a)  Tensor  storage.  The  lmth  component  of  W(l.m,y,y)  tensor  is  displayed  on  the 

electronically  addressed  spatial  light  modulator  (ESLM).  Lense  LI  images  W(l,m,x,y) 
onto  the  random  phase  function  exp  [;'0(.v,y)].  This  random  phase  function  is  provided 
by  a  high  resolution  E-Beam  fabricated  random  phase  code  transparency  (RPC).  The 
resulting  product  is  Fourier  transformed  via  lens  L2  onto  the  photorefractive  crystal 
while  a  reference  plane  wave,  deflected  by  an  angle  la  and  ma,  arrives  from  the 
acousto-optic  device.  The  superposition  of  the  first  beam  and  the  shifted  reference 
beam  creates  the  appropriate  interference  hologram.  The  space  bandwidth  product  of 
the  photorefractive  crystal  is  high  (109)  and  will  allow  the  storage  of  a  large  W(l,m,x.y) 
tensor,  (b)  Random  Phase  Coded  Matrix  Tensor  Multiplication:  When  G(x,y)  is 
presented  during  readout,  the  lm  correlations  of  W(l,m,x,y)  with  G(x,y)  sharpened  by 
the  random  phase  function  are  transmitted  along  the  various  (la,ma)  directions.  This 
beam  is  then  Fourier  transformed  by  lense  L3  and  captured  by  the  CCD  camera. 

After  the  initial  inputs  to  W(l,m.x,y)  are  made  the  two  phases  of  operation  will 
alternate.  An  image  will  be  input  and  the  correlations  determined.  Then  the 
appropriate  lmth  component  of  W  will  be  updated  by  an  amount  determined  by  the 
algorithm.  The  alternation  of  these  two  phases  provides  the  system  with  the  ability  to 
adjust  the  weighting  factors  and  to  learn  new  categorizations  as  they  appear. 
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Fig.  11a.  I-V  characteristics  of  a  typical  NMOS  pass  transistor  used  in  the  ESLM  circuit. 
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Fig.  1  lb.  Response  time  characteristics  of  ESLM.  The  upper  curve  is  plotted  on  a  longer  time 
scale,  the  lower  curve  shows  the  rise  and  fall  times  in  detail  on  a  shorter  time  scale. 


